Waves in Open Systems: Eigenfunction Expansions 
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An open system is not conservative because energy can escape to the outside. An open system 
by itself is thus not conservative. As a result, the time-evolution operator is not hermitian in the 
usual sense and the eigenfunctions (factorized solutions in space and time) are no longer normal 
modes but quasinormal modes (QNMs) whose frequencies u are complex. QNM analysis has been a 
powerful tool for investigating open systems. Previous studies have been mostly system specific, and 
use a few QNMs to provide approximate descriptions. Here we review recent developments which 
aim at a unifying treatment. The formulation leads to a mathematical structure in close analogy to 
that in conservative, hermitian systems. Many of the mathematical tools for the latter can hence be 
transcribed. Emphasis is placed on those cases in which the QNMs form a complete set for outgoing 
wavefunctions, so that in principle all the QNMs together give an exact description of the dynamics. 
Applications to optics in microspheres and to gravitational waves from black holes are reviewed, 
and directions for further development are outlined. 



I. INTRODUCTION 



A. Physical motivation 



Many concepts in physics rely on the idea of the normal 
pmodoQ (NMq) of a oonQorvativc ayotom. Thua one opoaka 



precisely the point of this paper: the familiar formalism, 
with very little alteration, apply to a wide class of these 
systems. 

The problem of EM waves escaping from an open cav- 
ity is brought into sharp focu s by advances in cavity 



of energy oigonotatoD, molecular orbitalo, energy bando, 
transitions between states and excitation energies (at the 
first quantized level), or of propagators in a Feynman di- 
agram (at the second quantized level). These concepts 
[for quantum and interacting (i.e., nonlinear) ayatcma arc 



198^ ; ^ 



QED: in Fabry-Pcr ot cavities ( Heinzen et al. 
Martini et at , |1987 ), i n superconductin g micr 
ities ( |Kleppnei) 7 l981 ; Goy et at , 198S), in s emiconduc- 



tor heterostructures (Yokoyama et al., 1990) and in mi 



crospheres which confine glancing rays by total internal 
reflectiori (see, e.g., Chang and Campillci| , 1996; Bandogh- 



rooted in eigenfunction Ojcpanoiono in tormo of the NMo 
of classical and Zmear systems; e.g., the photon propaga- 
tor in QED depends on the free classical electromagnetic 
(EM) field being expandable in plane waves. Such ex- 
pansions are possible because conservative systems are 
l aooooiated with hermitian oporatoro. 



dar et al, 19961 ). First, the spectrum is dominated by the 



resonances. 



Figure lb shows an experimental spectrum 
observed from a microsphere, similar in essence to Figure 
la. Second, the decay rate of an excited atom or molecule 



On the other hand, if energy can oooape to the out — 
side, the system would be open and nonconservative, the 
associated mathematical operators are not hermitian in 
the usual sense. The eigenfunctions are then quasinormal 
modes (QNMs) with complex frequencies ui. Figure la 
shows schematically the EM spectrum observed outside 
a linear optical cavity of length a, due to emission by a 
broad-band source, e.g., fluorescent dye molecules. The 
resonances, spaced by Alo k, ixcj a (where c is the speed 
nfliRbt), arc QNMs, which are rhpiractpristic of the ca 



1987 



De Martini et 
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Kleppner 




1981 




Heinzen ei 
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Barnes et. 


a, 1996 


)- 



the atom or the molecule of dimension 0.1 nm "knows" 
about its environment, on a scale of ~ 1 /im, even before 
any photon is emitted. 

Gravitational waves can be produced by matter falling 
into a black hole, and may soon be observed (see, e.g., 
199^ ). 



Abramovici et al. 



"a'^nty 



The curvature of the inter- 
vening space forms an effective potential akin to a cav- 
ity, from which the gravitatio nal waves eventuall y es- 
cape. Numerical simulation s ( Vishveshwara j_ 1970 ;_ De- 



rather than of the source The resnnance width -v is de- 



termined (in the a.bsence of ahsorption) by the amount of 
leakage. In the limit of zero leakage, these QNMs reduce 
to the NMs of a conservative cavity. It would be both in- 
teresting and useful if the physics of open systems could 
be discussed in terms of the discrete QNMs, in other 
words, an eigenfunction expansion for the dynamcis. 

This task is nontrivial because there are very few 
known results for eigenfunction expansion in nonconser- 
vative, non-hermitian systems. Yet the formulas we show 
in this article will look extremely familiar — and that is 



tweiler and Szedenits, 197£ ; 3mari , 1979; Stark and PI 



ran| , 1985| ; [Anninos etal. , 1993) show that the amplitude 
is dominated, at intermediate times, by a ringing signal 
(Figure 2) ^a^e'^'^^*, where j labels the QNMs, and 
each complex frequency ujj is characteristic of the back- 
ground geometry rather than of the emitting source. 

In these and many other examples, the resonances 
or QNMs often dominate. But how complete are the 
QNMs? In mathematical terms, this means, first of all, 
whether any function (pix) can be expanded in terms of 
the QNMs fj{x): 



1 
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<j>{x) 



3 



ajfj{x) 



(1) 



More importantly, do the resonances represent the dy- 



nan ics exac 



(2) 



for alH > 0? These expansion should make no reference 
to the environment or external bath. One would wish 
to establish conditions under which these expansions are 
valid, and in circumstances where they are not, to char- 
acterize the "remainder" . 

A projection formula, i.e., some sort of inner product, 
is needed to obtain the coefhcients aj from (j){x). Then, 
initial value problems become formally trivial: take the 
initial data, project out aj, and evolve by aj i— > aje"*"^*. 
Moreover, to the extent that similarities can be estab- 
lished with the conservative case, one should be able to 
transcribe the tools of mathematical physics, and to seek 
a parallel formalism — including, e.g., second quantiza- 
tion. These developments should lead to a unifying view 
of QNMs in different systems. 

This article discusses recent progress towards these 
goals. 

Open systems are a special class of dissipative systems, 
which have been stu d ied b y many au t hors, i ncluding 
Fcynman and Vernon (1963), UUersma ( 1966| ), paldeira 
and Lcggctt ( 1983 ). In all cases, the bath degrees of 
freedom are first included, and then eliminated from the 
path integral or equations of motion. 

Similarly, an open cavity (linear dimension a) can be 
embedded in a "universe" (linear dimension A —f oo), 
and the totality is a conservative system with the modes 
of the universe (MOU) forming a continuum {Alu oc 
A"-'^ ^0). A rigorous theory of lasi ng in 1-d cavities 
has be en developed using the MOU ( Lang, Scully and 
LanQ [l973a ; Lang and Scully , 1973b ) and the ideas can 
be generalized to higher dimensions, e.g., o ptics in micro- 
sph eres (see, e. g., Ching, Lai and Young , 1987a ; Ching 
Lai ind Young, 1987b), including nonlinear optical pro- 
cesses such as spontaneous and stimulated Brillouin 
scatteri ng ([Ching, Leung and Young , 1990; Leung and 
You0 |1991b| ). 

Since each eigenfunction or QNM is a resonance, the 
expansions (|l|) and (^) in effect expresses the idea of reso- 
nance domination. In cavity QED, re sonance domination 
was first discussed by Purcell ( 1946 ). He proposed that 
in the Fermi golden rule, the the density of states per 
unit volume po{oj) — uj'^/{tt^c^) should be replaced by 
p{oj) ~ D/{2'yV) for a D-fold degenerate QNM of width 
7 in a cavity of volume V. This leads to an enhancement 



factor oi K = p/po ~ {l/87r)DQ{\^ /V) for spontaneous 
emission, where A is the wavelength of light emitted and 
Q is the quality factor of the cavity. Purcell' s argument 
can be made more precise using the MOU ( Ching, Lai 



and Young, 1987a; Ching, Lai and Young, 1987b): the 
photon modes are re-distributed — enhanced at the res- 
onance and depleted away from the resonances. 

These works, in company with much earlier works 



in metastable states in nuclear physics (Gamow, 1928 
biegert , 1939| ) typically use some resonances for an ap 
proximate description. Our focus will be on the use of 
all the discrete resonances to obtain, if possible, an exact 
description — and where this fails, to characterize the 
remainder. 



B. The mathematical problem 

Waves in open systems may be governed by the 
Schrodinger equation (e.g., nuclear physics), the wave 
equation (e.g., optics) or the Klein-Gordon (KG) equa- 
tion with a potential V (e.g., each ang ular momentum 



secto r of linearized gravitational waves ( Chandrasekhar 



1991)). The wave equation can be mapped exactly into 



the KG equation, and as far as time-independent poblems 
are concerned, the KG equation can be mapped into the 



Schrodinger equation by relabeling lu 



Thus, many 



properties are similar (although we shall also highlight 
the key differences) , and we shall focus on the scalar wave 
equation, on a certain domain TZ, 



D $(r, t) EE [p{r)df - V^] $(r, t) = 



(3) 



This describes the scalar model of EM {p = dielectric con- 
stant, and henceforth c = 1) or elastic vibrations (p = 
density). In this article, we shall for the most part con- 
sider the 1-d version (V^ i-^ 9^), restricted to a half line 
< a; < oo; the "cavity" is the interval TZ = [0,a]. A 
node is imposed at the origin (a totally reflecting mir- 
ror), and waves escape through the point x = a {a par- 
tially transmitting mirror) to the rest of the "uni verse" 



m a, oo 



This model describes a laser cavity (Lang 



Scully and Laml , 1973a ) , as well as gravitat i onal r adi- 
ation from a stellar object (Price and Husain, 1992), or 



the transverse vib r ation s of a string ( Dckker| , 1985 ; 
Leung and Young, 1987). 



Lai 



For conservative systems, one usually imposes $ = 
on the boundary dTZ (2; = a in the 1-d case) ^. The 
eigenfunctions or NMs are factorized solutions $(r,i) = 
/(r)e~*'^*. The nodal condition on dTZ implies that TZ 
is closed, so that energy cannot escape. Mathematically, 
the operator — is hermitian; thus uj is real and the 
eigenfunctions {/} form a complete orthogonal set. 



^The Dirichlet condition can be replaced by the Neumann 
condition. 
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In contrast, if <I> satisfies the outgoing wave condition 
on dTZ; this defines an open system. The eigenfunc- 
tions, labeled by an index j, satisfy (in 1-d) d'^fj{x) = 
—ujjp{x)fj{x), with Im LUj < 0. Note a doubling of 
modes compared with the conservative case; If LUj is 
an eigenfrequency, then so is uj-j = — ^j, but since 
oj'^j ^ oj'j, f-j(x) and fj{x) are linearly independent. 
The question then is whether these {/j} are complete in 
the sense of (0) and (0). 



II. COMPLETENESS 

Even though QNMs do not always capture all the 
physics, it is still useful to first establish the baseline 
case in which the QNMs are complete on 7?. = [0, a]. This 
holds provided two conditions are satisfied: (a) p(x) has 
a step or stronger discontinuity at x = a (discontinuity 
condition); (b) p{x) = 1 (or other constant value) for 
X > a ("no-tail" condition). These two conditions can 
be understood heuristically. A discrete representation 
of (l){x) such as (|^) could at best be valid over a finite 
interval, for which a discontinuity provides a natural de- 
marcation. Moreover, a dynamical expansion such as (|^) 
would require that waves are not affected by the outside 
(a, oo); this can only be the case if p is trivial outside, 
and thus outgoing waves are not scattered back. 

We now sketch the main elements of the proof of 
completeness, which will also provide the framework for 
understanding other contributions when the QNMs are 
not complete. The Green's function G{x,y;t), defined 
by DG{x,y;t) = S{x - y)S{t), with G{x,y;t) = for 
< < 0. Its fourier transform G{x,y;uj) is given explic- 
itly as G{x,y,uj) = f{uj,x)g{u},y)/W{uj) ior Q < x < y. 
Here / and g are homogeneous solutions to the time- 
independent wave equation at frequency with / satis- 
fying the left nodal boundary condition (/(w, 0) = 0) and 
g satisfying the right outgoing wave boundary condition 
{g{io, x) oc e*"^ for x oo). Their Wronskian is W, and 
the combination fg/W is independent of normalization 
convention. 

For t >0, G{x, y; t) is evaluated by the inverse Fourier 
transform with the contour of integration closed by a 
large semicircle in the lower half w-plane. In general, 
there are three different contributions. 
Large semicircle — prompt response 

The contribution along the large semicircle {\lu\ — > oo) 
gives the short time or prompt response. For large t, the 
factor e~*'^* provides sufficient damping in the lower half 
cj-plane to control the asymptotic behav i or an d it can 



be shown (Bachelot and Motct-Bachelot , 1992; Leung ^ 
Liu and Young , 1994aD that the prompt response always 



vanishes for t > tp(x,y), where tp(x,y) is the geometric 
optics transit time between x and y. 

If the discontinuity condition is satisfied, one has a 
stronger result: by a WKB estimate, this contribution 
vanishes for all t > 0. It is natural that the \uj\ oo be- 



havior should be sensitive to short spatial length scales, 
in particular a discontinuity. 

Singularities in g — late-time tail 

The function f{uj,x) is obtained by integrating the 
time-independent wave equation, in which uj appears an- 
alytically, through a finite distance from to x. Crudely 
speaking, this is a finite combination of analytic func- 
tions of cj, so f(uj,x) is guaranteed to be analytic in uj 
(NewtoD, 19601 ). On the other hand, g{uj,x) could have 
discontinuities in uj (typically a cut on the negative Im uj 
axis), because it is obtained by integrating through an 
infinite interval starting from a; = oo. 

However, if the "no tail" condition is satisfied, then 
one can impose the right boundary condition for g{ijj,x) 
at X = a+ ; one then again integrates only through a finite 
distance, thus guaranteeing the analyticity of g as well, 
and removing this contribution. 
Zeros of W — QNM contributions 

Finally, there are contributions from the zeros of W{uj) . 
At a zero Uj of VF, the functions / and g are linearly de- 
pendent: fj{x) = f{ujj,x) — Cjg{ujj,x); thus fj satisfies 
both the left nodal and the right outgoing- wave boundary 
conditions, and is an eigenfunction or QNM with eigen- 
frequency UJj . 

In general, all three contributions are present, and cor- 
respon d nicely with the main features o f the ti me domain 



signal ( Ching et al. , 1995a ; Ching et al. , 1995c ) , as shown 



in the spacetime diagram in Figure 3. The contribution 
from the semicircle in the w-plane leads to the initial 
transients — a wave propagating directly from the source 
point y to the observation point x (rays la, Ih). The ze- 
ros of W give rise to the QNM ringing at intermediate 
times, which are caused by repeated scattering from the 
potential (ray 2). The cut in the w-plane, especially its 
tip near = 0, then gives the late-time tail: waves prop- 
agate from the source point j/ to a distant point x' , are 
scattered by p{x'), and return to the observation point x 
(ray 3). 

We now focus on the case where the discontinuity and 
"no-tail" conditions are satisfied. At each zero LOj, the 
residue is related to dW{ujj)/duj, which can be evaluated 
using the defining equations for / and g to be 



1 dW{LUj) 
Gj duj 



— ^ujj / p{x)fj{x) dx + ifj{a ) 





(4) 



Using the normalization (and phase) convention 
ifjlfj) — 2(jjj, one then obtains 



G{x,y;t) = 1^ l/^.(x)/,(y)e-^"^* (5) 



The initial conditions G{x,y;t = 0) = 0, 



February 7, 2008/ Waves in open systems: Eigenfunction expansion / Draft 



4 



p{x )f)tG{x, y;t = 0) = 6{x — y) lead to^ 



et ai, 1990| ), and in fact applies also to systems without 



Z, — iu^ 

3 •' 

^^^fA^)fo{y) = Kx-v) 



(6) 
(7) 



for x,y € TZ. The identity (|^) has an obvious analog in 
conservative systems; the factor 1/2 accounts for the dou- 
bling of modes. On the other hand (^) becomes vacuous 
in the limit e — > 0. 

The heart of the eigenfunction expansion is contained 
in (^), and The identity ^) obviously leads to 
(|I|), while (|) wiU be seen to lead to (|). 



discontinuities. 

To illustrate, let p{x) = 1 + + /i8(.T — b) for x < a 
and p{x) = 1 for x > a (i.e. a "dielectric rod" with a 
"bump"). Figure 4 shows the QNM positions in the em- 
plane for — 4.0, /i = 0.4, compared with /i = 0. The 
two are quite different, in that for the former case Im ujj 
are oscillatory in the mode number j. The perturbative 
results based on (||) and are also shown in Figure 4; 
although these formulas mimic the case of real frequen- 
cies for conservative systems, they nevertheless give the 
correct shifts, even for the imaginary parts. The sum over 
intermediate states is discrete, at a spacing Aw ~ tt/o; 
thus the effective expansion parameter is Awj ^ fia/n, 
a perspective not apparent in the MOU approach. 



III. TIME-INDEPENDENT PROBLEMS 

This simplest and most useful applications concern 
time-independent problems, especially time- independent 
perturbation theory. 

A. Perturbation 

Let p{x) — po{x) + p,V{x), where \p,\ ^ 1, and V{x) 
is nonzero only in the cavity. The system with po{x) is 
exactly solvable and its QNMs form a complete set for 
X £ [0, a] in the sense of (||), (||) and (0). Thus, it is 
straightforward to develop a discrete representation for 
the frequency shifts loj 



2 (2) 



(1) 



,(0) 



(8) 



,(2) 



(0) 



UJ, 



(0),^(0) 



UJ. 



(0) 



} (9) 



etc., where the matrix elements are Vji = /j Vfl dx 
These complex formulas give the shifts in the resonance 
positions and the widths — thus they contain twice 
the information compared with their superficially simi- 
lar conservative counterparts . The wavefunctions can be 
calculated in a parallel way ( Leung et al. , 19941: ) . 

The first-order correction has been known for a long 
time for the Schrodinger equation, but here the normal- 
ization condition (see (^)) is stated through a s urfac e 
term rather than via regularization (Zcldovich, 1960), 
which is less convenient computationally. The first- 



order result has been generalized to the EM case (Lai 



B. Physical examples 

Some interesting applications relate to optics; the in- 
terface between two media naturally gives a discontinu- 
ity, and vacuum outside the system naturally satisfies the 
"no tail" condition. Consider a microsphere of radius a 
and refractive index n. Glancing rays suffer total internal 
reflection and are confined, but evanescent waves cause 
some leakage, and this makes the system an open one. In 
practice, experiments are often done on droplets falling 
in air, slightly oblate (ellipticity e ^ 10~^ — 10~^) due to 
viscous drag. The breaking of spherical symmetry lifts 
the degeneracy among the 21 + 1 members of a multiplet. 
This proble m has been treated by brute force numerical 
calculation (Barber and Hill, 198S), but the perturbative 



formalism allows an analytic expression ( Lai et al , 199C) 



UJm 



3m2 



1(1 + 1) 



(10) 



in which m is the azimuthal quantum number, and 
e = [rp — re) /a, where rp and r^. are respectively the po- 
lar and equatorial radii, and a = {rpr"^)^^^. This pertur- 
bative formula has been used to interpret spectroscopic 
measurements and to determine the ellipticity e, either 
by measuring t he splitting between two neigboring lines 



(to and TO-f 1) (Chen et ai, 199l| ), or by determining the 
t otal spread of the multiplet (from to = to m = ±1) 



( phen et al , 1993). Two other developments along this 



line should be mentioned. 

In one interesting experiment 



(Arnold, Spock and 



Folan, 199C), a quadrupole field Eq imposed on a lev- 



itated droplet caused a spheroidal distortion with an el- 
lipticity e cx Eq/s, where s is the surface tension. The 



^Subject to the validity of term- by-term differentiation and 
the limit t ^ 0. 
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splitting of the multiplet was determined spectroscopi- 
cally and the ellipticity e was then found from (10), al- 
lowing s to be determined, down to length scales of tens 
of fiTa. (Hitherto, surface tension is known only macro- 
scopically, at length scales of say 1 mm or more.) 



simultaneous expansion of a pair of functions 



In another experiment (Swindal et al, 1992), light is 



launched into a slightly oblate microdroplet, propagating 
along great circles inclined at an angle 9 to the equato- 
rial plane. In photon language, the angular momentum L 
satisfies m/l = Lz/L = cosd, and precesses at the group 
angular velocity = dujm/dm = Lij\e\7n/[l{l + 1)]. The 
precession is observed as variations in intensity when a 
particular spot at an angle 9 is observed. The observed 
period of typically 50 ps then allows the ellipticity to be 
determined. 

On a very different length scale, the perturbation re- 
sults have also been applied to gravitational waves. The 
QNM spectrum of linearized gravitational waves propa- 
gating on a Kerr or Schwarzschil d background has been 
extensively studi ed analytically ( Chandrasekhar j_ 1991 ) 
and nu meric a lly ( Vishvcshwara , 1970 ; petweiler and Szc- 



deni^, |1979| ; ^mari] , |1979| ; [Stark and Pira4 |1985|), and 



their complex frequencies are known in detail ( Leaver 
1985|; [Guinn et al|, 1990| ; Leavci , 1992; Andersson and 



Lim ,aeu4 |1992 f ^ollert| , |1993 ). However, black holes are 
perturbed by interactions with its surroundings (e.g., a 
massive accretion disk). The shifts may provide a probe 
of the intervening spacetime curvature. This situation is 
obviously one to which perturbation theory can be ap- 
plied. 

The first-order perturbation theory is expressed by the 
KG analog to (P), further generalized to allow for a "tail" 



( [Leung et al\ , [1997] ). We here present some results for a 
Schwarzsch ild black hole pertu rbed by a thin static shell 
of matter ( Leung et al. , 1997). Figure 5 shows, in the 
complex w-plane, the locus of the lowest QNM for Z = 1 
scalar waves propagating on a Schwarzschild background, 
as a light shell (0.01 times the mass of the hole) is placed 
at different positions Vg. The first-order perturbative re- 
sults (dashed lines) capture the qualitative features of the 
exact results (solid lines). 



IV. DYNAMICS 



T ime-dependent problems require two independent 



3 



-iujjp{x) 



(11) 



using the same coefficients Oj 's for both components. The 
factor p{x) turns the second component into the conju- 
eate momentum. The use of two components is natura l 



(jFeshbach and Villars[ , [1957] ; |Unruli[ , [1974[ ; [For4 [1975[) 
but here the outgoing wave condition ^ (j){x = a+) = 



-0'(x 



") links them in a novel way. 



Compared with the case of NMs, twice the degrees of 
freedom {a.j and a_j) are used to satisfy twice the con- 
ditions {(f> and (j)). In the conservative limit, uP^^ = w| 
and f-j{x) = fj{x). Thus, upon grouping pairs of terms, 

H^) = Ej>o^j/i(2^)' ^i^) = Pi^)Y.o>QCjfj{x), where 
hj = Qj + a-j, and Cj — —iujj{aj — a-j). So in this limit, 
the simultaneous expansion of two functions (</>, </)) using 
the full set of NMs j — ±1,±2, • • • is equivalent to the 
separate and independent expansion of cj) and (j)/p using 
only the NMs j = 1,2, • • •. 

The solution to the initial value problem is then 

$(x,t)= / G{x,y;t)^{y) + dtG{x,y;t)p{y)(j){y) dy 
Ja L ^ 

(12) 

From (^2|), one can then obtain (|ll] ) and also (||) with aj 
given by 



fjiymy) + fjiy)m\ dy 

-/,(a)</>(a)l (13) 



sets of initial data: 



^(x,t = 0) and 



2w,- 



All reference to the outside of the cavity is removed, be- 
cause the integral on (a, oo) can be collapsed to the sur- 
face term, by making use of the outgoing condition on 
the initial data and the retarded condition on G. 

The elimination of the outside is the crucial step in 
obtaining a self-contained description of the cavity. 

The projection formula (|3|) suggests the definition 
of a generalized inner product, and with it a linear 
space structure analogous to that for conservative sys- 
tems (Leung, Tong and Young , 1997a; Leung, Tong and 



pix)dt^{x,t = 0). 



One is thus led to consider the 



Young , 1997b). Consider the linear space of outgo- 
ing function pairs {(j>,4'), denoted as ket vectors \<p) 



^With a "tail", the QNMs are not complete, but it is not 
surprising that this does not matter for the first-order shift. 

^These are assumed to vanish at infinity, representing waves 
emitted at a finite time in the past. 

^The last condition is specified at a"*" because the system 
may have discontinuities at a; = a. 
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(0(a;),0(x))T. For a QNM, fj{x) 



-iujjp{x)fj{x), function p{x,u>) mu st satisfy th e usual Kramers-Kronig 

|l994b| ). 



where ujj is the eigenvalue. Time-evolution can be writ- dispersion relation (Leung et al. 
ten in Schrodinger form dt\(f>) = —iTi.\(p), where 



n 



p{x)- 







(14) 



The first component of the Schrodinger equation repro- 
duces the identification of cj) as p{x)dt^. The eigenfunc- 
tions, or QNMs, can now be defined simply by 7i — 

The projection (|l^) suggests a generalized inner prod- 
uct between two vectors \cf>) and \ip)- 



{il,\(t>) 



{ip(f> + %lj(f>)dx + 'ijj{a)(l)(a) 



(15) 



which is symmetric and linear in both the bra and ket 
vectors (rather than conjugate linear in the bra vector). 
The inner product of a QNM with itself agrees exactly 
with the generalized norm (|). Mor cover, the coefficients 
( |l3|) for the eigenfunction expansion can now be written 
compactly as 



(16) 



Under the generalized inner product, the Hamiltonian 
is symmetric: (-01 {n\<j)) } = {<j)\ {n\ip) } ee {ip\n\<j)). 
In the proof of this statement, there is an integration by 
parts; surface terms are incurred because the functions 
do not vanish at the end-points, but these are compen- 
sated exactly by the surface terms in the definition of the 
inner product. 

There is thus an almost complete parallel with con- 
servative systems, and the mathematical structure is in 
place to carry over essentially all the familiar tools based 
on eigenfunction expansions. For example, the symmetry 
of TL (similar to the hermitian property in the conserva- 
tive case) immediately leads to the orthogonality of the 
QNMs and hence the uniqueness of the expansion (11). 
The only exception is the lack of a positive-definite norm, 
and with it a simple probability interpretation — hardly 
surprising since probability (or energy) is not conserved 
in TZ. For example, the perturbation theory can now be 
obtained by an almost trivial transcription of textbook 
derivations. 

The formalism generalizes to a full line (— oo < x < 
oo); the expansion is valid in an interval [01,02], pro- 
vided (a) there are at least two discontinuities, with the 
extreme ones at oi and 02; (b) there is "no tail" on either 
side, and (c) the outgoing condition is imposed at the two 
ends o^ and o^ . Much the same formalism applies to the 
KG equation. Discontinuities must be present in the po- 
tential V{x), and the "no tail" condit i on is that V { x) — 
on the "outside" ( |Ching et alj |l995b| ; |Ching et a/.| , |1996D . 
Absorption can also be included; e.g., in the scalar model 
of EM, the complex and frequency-dependent dielectric 



To apply to cavity QED, it is necessary to (a) go from 
free fields to interacting fields, and (b) go from classical 
fields to quantized fields. Here we indicate the rudiments 
of this program. 

Consider a nonlinear wave equation: + \{x)^^ = 
where A(x) has support only within TZ (interaction oc- 
curs only inside the cavity). Expanding this in terms 
of the eigenfunctions of noninteracting theory, {d/dt + 



lOJi Oo 



-i /{2uj j)Y,kim'^3kimakaiam, where Xjkim 



Iq ^fjfkfifmdx Time-dependent perturbation of the 
linear theory (i.e., adding a time-dependent piece to p{x)) 
can be handled in very much the same way. The point 
is that the dynamics is now expressed through discrete 
variables aj{t) 



Since (/)(x) and (l){x) satisfy 



4'{x),(t>{y) 



iS{x — y) at 
commutation relations 



The commutators differ slightly from ujj , 



equal times, one obtains from ( 
for [oj , o^ 

the difference being a measure of the dissipation. Feyn- 
man rules can be developed in the standard way. This 
description, unlike the MOU approach, will be particu- 
larly useful when the cavity is only slightly leaky and 
the connection with a totally enclosed cavity is to be em- 
phasized. Applications to emission, absorption and other 
quantum phenomena are bei ng deve l oped. An important 
result will be to put idea of PurccU ( 1946 ) on resonance 
enhancement of decay rates on a firm footing, by writing 
it as the contribution of one discrete QNM, in a precise 
manner. 



V. NON-RESONANT CONTRIBUTIONS 

We have so far concentrated on the baseline case in 
which the discontinuity and "no tail" conditions are sat- 
isfied, and the QNMs are complete. However, in other 
cases there are non-resonant contributions, and a brief 
discussion is given in this Section. 

Gravitational radiation from Schwarzschild black holes 
is described by the KG equation; the potential V{x), re- 
lated to the background metric, is continuous and goes 
asymptotically as a;~^logx as x ^ -|-oo (apart from a 
centrifugal barrier) , violating both the discontinuity and 
the "no tail" condition. We indicate briefly how this sit- 



uatio n affect s the dynamics (Ching et 



et alj |l995cD . 



1995a 



Ching 



First, because there is no discontinuity, there is a 
prompt signal due to the contribution from the large 
semicircle. Secondly, the auxiliary function g(uj,x) has 
to be integrated from x — s- +00. This leads to a cut on 
the negative Im oj axis, controlled by the spatial asymp- 
totics of V{x). The cut extends all the way to w — 0, 
and determines the late time dynamics. Generically, if 
V{x) ~ x-°'{\ogxY, then (j){x,t) ~ t-^'{\ogtY, where p, 
and V are determined by a, /? and the angular momen- 



February 7, 2008/ Waves in open systems: Eigenfunction expansion / Draft 



7 



turn /. For example, for /3 = or 1, in general fj, = 21 + a 
and V = (3. It turns out, interestingly, that the case of 
the Schwarzschild black hole (a = 3, /5 = 1) is excep- 
tional for I > 0: the leading term ^^(2'+") \ogt vanishes, 
and the dynamics is dominated by the next leading term 
^-{2i+a) ^ Figure 2 shows precisely such a power-law tail. 
However, the signal at intermediate times would still be 
dominated by the sum over QNMs, so that in practice, 
a restricted and approximate notion of completeness still 
holds. 

For the Schrodinger equation, outgoing waves are de- 
fined by (f){x) ^ e^^^ as a; ^ oo, where k — \/2tJ instead 
of A: = w, and moreover fc > 0. Thus a continuum con- 
tribution cannot be avoided, and the system decays by 
a power-law at large times, due to the contribution 
ne ar threshold ( fc « ). This pheno mena was first noted 
by 
te 



^ iKhafinI (|19 
i) [|1961| 



5?? 



) and Winter (1961). In particular, Win- 
considered a 1-d well bounded by an impen- 
etrable wall and a delta-function potential and showed 
that the wavefunction inside the well decays as t"'^/^ as 
large t. However, if the potential is unbounded from be- 
low, no finite threshold exists, and the large t behavior is 
again given by the QNMs (Suen and Young, 1991). This 
turns out to apply to some models of mini-superspace 
( in which x is essentiall y the scale factor of the universe) 
( ^uen and Young , 1989), and there is the intriguing possi- 
bility that wavefunction of the universe at the end of the 
quantum era is given by the lowest QNM, independent 
of cosmological initial conditions. 

The power-law decay in the case of the Schrodinger 
equation comes from the threshold, whereas the power- 
law in the KG case comes from the spatial tail of V{x); 
the origins of the two are entirely different. 

Despite these complications, in time-independent 
problems, one is free to relabel uP' i~+ w, so the 
Schrodinger case should not be different. An applica- 
tion in this restricted time-independent sense is given in 
Leung and Yoimg (1991a). 



VI. CONCLUSION 

We have considered linear, classical waves in an open 
system of a certain type, and sketched a formalism in 
terms of an QNM expansion, in a manner analogous to 
the description of conservative systems by its NMs. Dis- 
sipation is contained in these discrete QNMs themselves, 
especially in Im Uj. Remarkably, apart from this fea- 
ture, almost nothing needs to be changed, and the fa- 
miliar tools for hermitian systems can be carried over. 
These results are nontrivial, in that if either the discon- 
tinuity condition or the "no tail" condition is violated, 
the QNMs would not be complete. The most significant 
feature in these cases is a power-law tail in the large t 
behavior. 

With these mathematical concepts in place, the way is 
now open for formulating QED phenomena in an open 



cavity, using a discrete basis. 
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FIGURE CAPTIONS 

Fig. 1 (a) Schematic diagram of the EM spectrum ob- 
served outside a hnear optical cavity of length a, 
due to emission by a broadband source inside the 
cavity, (b) Experimental fluorescent spectra (up- 
per) of a dye-doped cthanol microsphere. Reso- 
nance peaks are found to correlate well with the 
computed resonant mode positions (lower). Mode 
positions are plotted as half arrows in a horizontal 
plane and projected by broken lines. [Figure repro- 
duced with permission from Fig. 5 of Chapter 4 of 
Optical Processes in Microcavities ed. R. K. Chang 
and A. J. CampiUo, World Scientific (1996).] 

Fig. 2 Numerical simulation of the a mplitude of lin- 
earized gravitational waves propagating on a static, 
spherically symmetric blackhole background as a 
function of time. Note particularly the domination 
by a ringing signal due to the QNMs at intermedi- 
ate times. 

Fig. 3 A spacetime diagram illustrating heuristically 
the three different contributions to the Green's 

fimction. Rays la and lb show "direct" propa- 
gation without scattering, and correspond to the 



prompt response. Ray 2 suffers repeated scatter- 
ings at finite x' . For "correct" frequencies, such 
multiple reflections add coherently so that reten- 
tion of the wave is maximum. The amplitude of 
the wave decreases in a geometric progression with 
the number of scatterings, and thus decreases ex- 
ponentially in time. These are QNM contributions. 
For a wave from a source point y reaching a distant 
observer at x, there could be reflections at very 
large x' by V{x'), as indicated by ray 3. For a gen- 
eral class of V, such ray contributes to the late-time 
tail. 

Fig. 4 Positions of the QNMs for the dielectric rod 
model with a small perturbation (crosses) com- 
pared to the unperturbed ones (triangles). First- 
order (squares) and second-order (circles) pertur- 
bative results are also shown. 

Fig. 5 The locus of the lowest QNM for / = 1 scalar 
waves propagating on a Schwarzschild background, 
as the shell of magnitude /i = 0.01 is placed at 
different positions measured by r^. The solid line 
shows the exact results, and the dashed line shows 
the first-order perturbative results. The square cor- 
responds to a shell position of Vg/Ma = 2.22 (the 
extreme value that satisfies the dominant energy 
condition). The triangles (on the solid line) and the 
crosses (on the dashed line) refer to rg/M^ from 6 
to 60 at intervals of 6. The first-order result is seen 
to capture the qualitative features. 



9 



